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Abstract 

Freeman Dyson has questioned whether any conceivable experiment in the 
real universe can detect a single graviton. If not, is it meaningful to talk about 
gravitons as physical entities? We attempt to answer Dyson's question and find 
it is possible concoct an idealized thought experiment capable of detecting one 
graviton; however, when anything remotely resembling realistic physics is taken 
into account, detection becomes impossible, indicating that Dyson's conjecture 
is very likely true. We also point out several mistakes in the literature dealing 
with graviton detection and production. 

PACS: .03.65.Sq, .04.30.Db, .04.30.Nk, .04.80.Cc, .04.80.Nn, 95.30.Cq 
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1 Introduction 

The search for gravitational waves, one of the central predictions of general relativity, 
has been ongoing for several decades. Strong indirect evidence for their existence 
comes from the timing of the orbital decay rate of the binary pulsar PSR1913+16, 
and with the completion of LIGO II and other gravity-wave observatories, researchers 
expect a direct detection. Gravitons themselves present a thornier issue. Although 
physicists routinely speak as if bosons mediating the gravitational force exist, the 
extraordinary weakness of the gravitational interaction makes the detection of one 
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gravitational quantum a remote proposition. Recently, Dyson [Tj has suggested that 
detection of a single graviton may in fact be ruled out in the real universe. If so, is 
it meaningful to talk about gravitons as physical, or do they become metaphysical 
entities? 

In this note we attempt to answer Dyson's question "in principle." For both phys- 
ical and philosophical reasons the matter turns out to be not entirely trivial, and 
both considerations require that the rules of the game be defined at the outset. We 
concede at once that there appear to be no fundamental laws disallowing the detec- 
tion of a graviton, and so we take the approach of designing thought experiments 
that might be able to detect one. Such an "experimental" approach, however, also 
has ambiguities; allowing an infinite amount of time for an experiment in any imagin- 
able universe would certainly allow the detection of anything. Thus we impose a few 
physical restrictions: We consider only experiments using standard physics in four 
dimensions that could be performed in something like the age of the real universe. 

One must also agree on what sort of experimental result one would accept as 
establishing the existence of a quantum of gravity. When one recalls that to establish 
the photon picture of light took a decade, if not two, one sees that this "philosophical" 
issue requires perhaps even more negotiation than the physical one. We discuss the 
matter in greater detail in §|3]and £0 Moreover, in addressing Dyson's conjecture, one 
is inevitably led to rather subtle physics, which apart from its own intrinsic interest, 
has caused some confusion and errors in the literature. For pedagogical reasons, we 
attempt to explain the basis for our calculations as well as the others, leaving however 
most of the fine details to the Appendix and a more technical companion paper [2], 
henceforth BR. 

Should Dyson's question be answered in the affirmative, that is, if one decides that 
detection of a single graviton is physically impossible, this immediately raises the issue 
of whether it is necessary to quantize gravity. There are of course theoretical reasons 
for wanting to do so — after all, every other fundamental field is quantized — and there 
has been some discussion in the literature about the consistency of gravitational 
theories without quantization [3j HJ [6] . These discussions appear to be inconclusive 
and we do not enter into them. Rather, we restrict ourselves to the more limited 
question we have already set forth: can one experimentally detect a graviton? 

To begin, we assume the most general physics and units in which c = 1, but 
G = G and h = h, and write down the criterion to detect a graviton. 

2 Detection Criterion 

Regardless of the exact nature of the detector or interaction, we demand that to detect 
a single graviton with high probability, its path-length A in the detector should exceed 
a mean-free- path, or 

na\ > 1 , (2.1) 
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where a is the interaction cross section and n is the density of detector particles. We 
wish to maximize the left hand side of this expression. For simplicity we consider a 
detector of atomic hydrogen, which implies n = Md/(m p V), with Md the detector 
mass, V its volume and and m p the proton mass. Assume the detector is spherical 
with a diameter i = X max . Then, if a total number of gravitons iV 7 is incident over 
the lifetime of the experiment, the criterion for detecting a single particle becomes: 

-^iV 7 >l. (2.2) 

iV 7 can be expressed in terms of the source's graviton luminosity L 7 and distance R 

as 

Here, A d = ii£ 2 /6 is the effective cross-sectional area of the detector (i.e., the cross 
section averaged over possible paths), t s is the time the source is on and e 7 is the 
graviton energy. Thus, 

47r m p Hr e 7 

Reasonably, r s < M s /L by energy conservation, where M s is the mass of the source 
and L its total luminosity. We also expect that any astrophysical source emits only 
a small fraction of its energy as gravitons. Letting f 1 = L 1 /L yields 



/ 7 a M d M s 



— > 1 . (2.5) 



An m p R 2 e 

Thus far this is almost completely general. We now assume that the detector should 
not be so large as to undergo gravitational collapse to greater than atomic densities. 
It is well known [7J that the requirement that gravitational forces do not overwhelm 
the electrostatic forces supporting the object gives a maximum detector mass of about 
the mass of Jupiter: Md ~ (a/a g ) 3 ^ 2 m p , where a = e 2 /h and a g = Gm 2 /h are the 
fine structure and gravitational fine structure constants, respectively. And so, finally, 
the criterion for detecting a graviton becomes: 

The crux of the problem is to determine a for a given process and / 7 . 



3 Cross sections: Heuristic Arguments 

Whenever dealing with quantum gravity, it is fruitful to enlist analogies with ordi- 
nary quantum mechanics and electromagnetism, and we shall do so throughout. In 
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choosing an experiment, and hence a cross section to calculate, we also follow this 
stratagem and attempt to learn from history. A century ago two experiments were 
crucial in establishing the particle nature of light: the photoelectric effect and Comp- 
ton scattering of photons by electrons. The photoelectric effect of course presented 
severe difficulties for classical electromagnetic theory, which predicted that the en- 
ergy of the photoelectrons should increase with the intensity of the incident light, but 
should be independent of the light frequency, neither of which proved true. Einstein 
explained everything by introducing light quanta and postulating the celebrated one- 
liner: E = hu — W, which gave the photoelectrons' energy in terms of the frequency 
of the quantum and the work function, the energy needed for the electron to escape 
the surface of the detector. The number of photoelectrons simply became equal to the 
number of incident quanta. Whether the formula actually held, though, was debated 
for ten years as physicists plotted E versus v to get the slope h, but in 1916 Millikan 
was able to announce that "Planck's h has been photoelectrically determined with a 
precision of about 0.5% to have the value h = 6.57 x 10~ 27 " [S|. 

With perfect hindsight, a devil's advocate would point out that the photoelectric 
effect is not so conclusive. In semi-classical derivations of the photoelectric cross 
section such as Schiff's[9], the electromagnetic field is never quantized. Thus, 
although the computed cross section gives good agreement with experiment, nothing 
in the mathematics implies a photon. There are two answers to this rather subtle 
point: one is that the photoelectric effect is a fact independent of the derivation of 
the cross section, and it still can't be explained by classical physics. The other is 
that we will accept the existence of a quantum of light (or gravity) if the detector 
can "fire" when there are fewer than one quantum in the detector at a given instant. 
We make this statement more precise in £0 after we have estimated expected fluxes, 
but will consider a "gravito-electric" experiment in which gravitons knock electrons 
from a detector to be a viable candidate. 

Compton scattering might be considered a more conclusive demonstration of the 
existence of photons, because the change in wavelength under scattering is calculated 
by treating the X-rays as nothing more than billiard balls with quantum-mechanical 
energy E = hu and momentum p = h/X. Thus, we can also imagine a "gravito- 
Compton" experiment in which a graviton is scattered off an electron and one at- 
tempts to detect the electron's recoil energy. (Detecting the scattered graviton itself, 
in analogy to the usual detection of Compton-scattered photons, is not kosher because 
it only returns the original problem of detecting the graviton.) 

The gravito-Compton effect presents its own subtleties. The ordinary Compton 
scattering cross section is oc ~ (e 4 /m 2 ), where e and m are the electron's charge and 
mass. One might naively assume by the dimensional substitution e — > \fGm that 
gravito-Compton scattering should have cross section o ~ G 2 m 2 . Indeed, numerous 
authors in the 1960s and 1970s considered scattering of gravitons off neutral scalar 
particles. Dewitt, for example [lOj. gives for the nonrelativistic limit the differential 
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scattering cross section 

da G 2 m 2 



dQ sin A [ 



cos s (^) + sin\^ 



(3.1) 



with an analogous relativistic limit. Perhaps to one's surprise, although the G 2 m 2 
dutifully appears, (13.11) resembles the Rutherford scattering cross section with its 
9 = divergence rather than the Thomson cross section 8 ^-(e A /m 2 ), which is the 
nonrelativistic limit of the Compton cross section. This is a hint to beware analogies 
between electromagnetism and gravity. To the best of our knowledge the scattering 
of gravitons off true fermions has not been calculated, but conversion of gravitons 
to photons by scattering off charged scalar particles has [11J, with the result that for 
unpolarized incident radiation: 



da e l G 



dQ 4vr 



cot 2 



sin z 2<p + cos 2 9cos 2 2<p . (3.2) 



This expression also diverges, but is larger than (13. ip by a factor of a/a g (m p /m) 2 ~ 
4 x 10 42 . As to why these cross sections diverge, the usual statement is that it is 
due to the long-range nature of the gravitational and electromagnetic force, but the 
different forms of the divergence or lack thereof in ordinary Compton scattering show 
that more detailed considerations play a role. 

Because the total cross sections of (13.11) and (13.21) are formally infinite, in practice 
one must put in a cutoff in the 9 integration. This can be accomplished by relating 
9 to the initial graviton and final electron energies in the usual Compton scattering 
formula; however certainly for (13.11) any such detail is irrelevant because G 2 m 2 ~ 
lCT 110 cm 2 . In other words, G 2 m 2 ~ lO -45 ^, where i pi ~ lCT 33 cm is the Planck 
length. Regardless of what other numbers one inserts into (12.61) . the detection criterion 
cannot be satisfied. Gravitons cannot be detected by Compton scattering of gravitons 
off neutral particles. 

All hope is not lost, however. One might reasonably assume the enormously larger 
(13. 2p must be the nonrelativistic result for true electrons and contemplate using it 
as the basis for an experiment. In fact no speculation is required. For any "gravito- 
atomic process," the actual cross section turns out to be, within numerical factors, 
the square of the Planck length, in our units GTi ~ 10~ 66 cm 2 , comparable in size to 
(13. 2p but without divergences. To get some idea of how this comes about, consider an 
idealized gravitational wave detector, consisting of two balls on a spring that behave 
as damped harmonic oscillators under the passage of a gravitational wave: 

£ + £/r + uj 2 ^ = acceleration, (3.3) 

where £ is the displacement of the balls. The basic definition of differential cross 
section is: 

da = Pjergs/s) ^ 
I{ergs/ s/cm 2 ) 
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By computing via (I3.3P the power absorbed in the detector, P, and dividing it by 
the incoming flux I ~ h 2 uj 2 /G of the gravitational wave (see §|5]), is not difficult 
to show [12] that near resonance ± u \ « u„) the cross section is given by the 
Lorentzian 

= TrGML 2 (u 2 /T )sin*e 
° (M-^) 2 + (l/2r ) 2 ' { ' } 

The mean cross section is just J adu, which gives a avg ~ GML 2 u Q . If we now 
assume that L = a = h 2 /me 2 , the Bohr radius, M = m, the electron mass, and 
impose the Nicholson-Bohr quantization condition mua 2 = h, then we immediately 
have a avg ~ G%, the square of the Planck length. Notice that the "atomic process" 
enters only in determining the size of a; if the atom were bound by gravity, the 
Nicholson-Bohr condition would still apply, but but a would be about 10 39 times 
larger and u about 10 78 times smaller. 

Consequently, the "Planck-length squared" cross section stems directly from the 
imposition of angular momentum quantization onto the detector and nothing more. 
This result has the immediate consequence that it is physically impossible to detect 
a single given graviton with high probability. Detection criterion (12.11) can now be 
written approximately as 

for detector mass M and size I. However, M/t < 1 for any object that is not a 
black hole. If we take i v \ = 1, then m p = 10~ 19 . The smallest conceivable detector is 
presumably about the size of a proton, or i ~ 10 20 , making the inequality impossible 
to satisfy. Smolin [6] has argued that this conclusion holds for any physical detector 
and that therefore it is impossible for gravitational radiation to be in thermodynamic 
equilibrium with it surroundings. 

However, this still leaves open the possibility that given a large flux of gravitons, 
a suitable detector might be able to detect one or more of them. Therefore, heartened 
by a non-divergent cross section that is forty some orders of magnitude larger than it 
could have been, we proceed to more detailed considerations. 



4 Matrix Elements 

The near universal prescription for the quantum mechanical calculation of cross sec- 
tions is to select an interaction Hamiltonian between the particle and field, then 
employ perturbation theory to calculate the transition probability between indi- 
vidual states: 

Pab oc — | < b\H\a > \ 2 



h 
T 1 



$>tH* a d 3 x\ 2 , (4.1) 
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where \l/ a and are the initial and final wavef unctions. Multiplying by the density 
of final states p = dn/dE gives the famous golden rule and the total transition rate 

r = =^-\<a\H\b>\ 2 p. (4.2) 

The differential cross section da for the process under consideration then follows 
immediately from the fundamental definition (13. 4p as the ratio between TlujT and the 
incoming flux, I. We will follow this prescription, but for gravity it turns out that the 
determination of the Hamiltonian and consequently V does involve subtleties, and we 
will try to clarify them. Doing so not only proves necessary to get the correct a but 
presents a few surprises. 

Consider first the calculation of the matrix element for the electromagnetic field via 
the standard semi-classical prescription. In this approach one begins by writing down 
the classical interaction Hamiltonian H = — e/m A-p, where A is the vector potential 
and p is the electron momentum. To compute | < b\H\a > | textbooks discuss replac- 
ing p by the momentum operator —ihV, but in fact what is usually done is to apply 
the correspondence principle in reverse: —ihV — > p — > mdr/dt, by which point both 
h and quantization of the field have dropped out of the problem. Additionally, r in 
the "energy representation" is taken to have time dependence e iu}abt , where hu^ is the 
energy difference between two levels a and b, implying that p = imc^r. One more- 
over assumes that the field itself is varying sinusoidally, i.e., A = A e l( ^*~ k ' r < ) . The 
time-independent part of the matrix element thus becomes A euj a b\ < 6|e* kr r|a > |. 
Usually, the wavelength in the problem is much larger than atomic dimensions, allow- 
ing one to make the "dipole approximation," e* k r = 1. And so, the matrix element 
one actually computes is simply | < b\H\a > \ — A eu ab \ < b\r\a > \. 

For the gravitational case, we copy the procedure, considering the interaction of 
an electron with a gravitational wave. A plane gravitational wave can be viewed as a 
weak perturbation h^ v traveling on a flat background T]^, such that the full spacetime 
metric becomes: 

In analogy to the electromagnetic case, we take for a monochromatic gravitational 
plane wave = he^e 1 ^ 1 "^, where is the polarization tensor. Because the 
perturbations are assumed weak, one takes the amplitude h « 1. 

To illustrate how to get an interaction Hamiltonian, it is perhaps easiest to proceed 
as follows: Recall that in special relativity, the Minkowski metric can be written dr 2 = 
—ri iiV dx iL dx v = dt 2 (l—v 2 ). This leads to the free-particle Lagrangian L = —m\/l — v 2 . 
For an interaction with a gravitational wave we replace rj^ by = rj^ + to get 



L = —mJl — ViV % — h fiu v^v' / , 



(4.3) 
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Here, = dx^/dt; we use the summation convention throughout and Latin indices 
refer to spatial components. Although this Lagrangian is not manifestly covariant, it 
serves for our purposes. 

Taking dL/dv a yields by definition the canonical momentum n a . In the nonrela- 
tivistic limit v « 1. To lowest order in v the canonical momentum becomes: 

where 5 l a is the Kronecker delta and p a = mv a is the ordinary momentum. By 
definition, the Hamiltonian is H = 7r a v a — L. Working again to lowest order gives: 

H _ Pif_ 1 h a pp a p p 
2m 2 m 

(Here we have ignored a physically irrelevant constant m.) If we define T^ u = mv^v u 
to be the classical stress-energy tensor for an electron (integrated over volume) the 
interaction Hamiltonian is then 

H int = \h^ v , (4.4) 

and we see that T^ u plays the same role as p did in the in the electromagnetic case, 
while h^ u corresponds to A. This derivation agrees with the more general result found 

in [la nil. 

Now, the polarization tensor e^ v has only two independent components. For the 
moment let us make the standard choice en = — e<i2 and eyi = e%\. If we again make 
the dipole approximation, then sandwiching H int between the initial and final state 
functions gives the required matrix element 1/2 < b\e^PiPj\a >. To evaluate this ex- 
pression exactly, one should let p — > —ihd/dx. However, if following the semi-classical 
procedure, we set p = mx, then we have m 2 < b\e l ^XiXj\a >, which dimensionally is 
m 2 uj 2 < b\e l ^XiXj\a >. Indeed, for the harmonic oscillator wavefunction, the expecta- 
tion value m 2 u 2 < a\x 2 \a > is identically equal to < a|p 2 |a >, and one can also show 
that the matrix element m 2 < b\xiXj\a >= m 2 uj 2 < b\xiXj\a > for any case in which 
all the transition frequencies u are the same. Thus we expect: 

<b\H\a> w ^^J^te ij x i x j ^ a d 3 x 

V ( 4 - 5 ) 

mhuj ij n 
— 2 % 3 ' 

where 

Dij = J ^bXiX^ad^X. (4.6) 

In principle, all we need to do now to get V is insert the relevant wavefunctions 
into Dij and calculate | < b\H\a > \ 2 . The matter is not so simple, however. First, for 
any absorption or scattering process the graviton can be incident from any direction. 
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Therefore, we need to square (14.51) for each polarization, average the results over the 
unit sphere and sum them. To do this, take |ey| = l/v2. Then, 



< | Ail 2 >= — / (19. 

47T 



2\m b xy^ a \ 2 + \\^ b {x 2 -y 2 )^ a \ 2 



< |Ail 2 >= -(AjAj - ttI^I 2 )- (4-7) 



2 

In BR, we give a quick way to evaluate this expression. The result is: 

2 1 

-(AjAj - - 

In what follows, we will thus assume that 

| < b\H\a > | 2 = < \ A/ >, (4.8) 

where < |Aj| 2 > is given by Eq. (14.71) . The extra 1/4 in this equation gives an exact 
expression in a sense explained in the Appendix. It arises, in part, when demanding 
that detailed balance be satisfied among the absorption, spontaneous and stimulated 
emission rates. The appearance of this factor also helps resolve some discrepancies in 
the literature and leads to some fundamental issues regarding gauge transformations 
in general relativity. Because these considerations are somewhat more technical than 
the general discussion and do not seriously affect the numerical result, we leave greater 
detail to the Appendix and to BR. Here, we proceed to calculate a cross section. 



5 Ionization cross section 

Despite our devil's advocate insistence that nothing in the previous discussion implied 
the existence of quanta, we again defer comment until §[7] and choose to calculate 
the gravito-electric (ionization) cross section of atomic hydrogen. Although some 
metals have work functions that are somewhat smaller than hydrogen's, hydrogen's is 
typical and the element does compose most of the universe; for this reason we chose 
a hydrogen detector in §2J It turns out that in the energy regime of greatest interest, 
the result is not sensitive to the work function. 

As in the standard calculation of the photoelectric cross section[9], we take the 
initial state to be the ground (Is) state of hydrogen and for the final state a plane 
wave with momentum k: 

*" = ^~ r/ ° ; *' = z^ ,k " r - (5 ' 1) 



where L is the normalization constant. The use of the plane wave final state (the 
Born approximation) is reasonably good for large final momentum, or a 2 k 2 >> 1, 
which is again the case of most interest. 
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Before one computes the matrix element with these wavefunctions, however, there 
is a complication. As in the ordinary photoelectric effect (see Schiff, p. 287), one 
should average the direction of the incident graviton over all angles. This average, 
though, is exactly the one we have already performed, leading to Eq. (14.71) . and 
therefore, in evaluating the matrix elements expression (14 .8p should be used for 
\<b\H\a>\ 2 . 

This actually simplifies matters greatly, because with the gravitational field aver- 
aged over all directions, we may assume that the electron is ejected along the z-axis 
and that 9, the angle between k and r, is the usual polar angle. Once again the 
integrals can all be carried out analytically in terms of elementary functions and the 
result is: 

I < b\H\n > I 2 — S^^frWaVfc 4 ) 

\<b\H\a>\ -— L3(1 + a 2 P) 8 W 

The total absorption rate T follows directly from the golden rule (14.21) . With the 
standard expression for the density of states p = L 3 /(2irh) 3 kmdQ, this gives 

_ 3(32) 2 h 2 uj*m 3 a u k 5 



4-5vr h 3 (l + a 2 k 2 ) 



; dn (5.3) 



The differential cross section now follows from (13 .4p as da = huT/I. For a gravita- 
tional plane wave, the equivalent of the Poynting flux is[12] I = a/Z^e^e*- 7 /16irG = 
uj 2 h 2 /8ttG, yielding 

6(32) 2 Gu 3 m 3 a ll k 5 ,„ 

da = — — - — ^ — ——dVl (5.4) 

5 Ti 2 {l + a 2 k 2 f V ; 

But the photoelectric equation requires E = (hk) 2 /2m = hu — e 2 /2a, where e 2 /2a 
is the hydrogen work function. Hence, 1 + k 2 a 2 = 2muja 2 /Ti. Inserting this into the 
above and integrating yields, finally, the total cross section 

3072 7T hG(ak) 5 
a ~ 5 (l + a 2 k 2 f (5 - 5) 

As promised, a ~ TiG to "within factors of order unity," but is strongly dependent on 
ak. (Indeed, for ak » 1, the dependence is the same as in the ordinary photoelectric 
effect.) We attempt to evaluate this factor in the following section. 



6 Graviton Production 

The second half of the task is to determine the graviton luminosity of astrophysical 
sources, that is estimate / 7 in Eq. (l2.6p . 

There are only a few conceivable sources of gravitons: spontaneous emission of 
gravitons from neutral hydrogen, black hole decay, bremsstrahlung from electron- 
electron collisions in stellar interiors and conversion of photons to gravitons by inter- 
stellar magnetic fields. We examine each in turn. 
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Spontaneous emission of gravitons, such as that from the 3d — > Is state of hydro- 
gen discussed in the Appendix, will not produce gravitons energetic enough to ionize 
hydrogen. One can, alternatively, imagine a gravito-Compton experiment with the 
cross section o ~ e 2 G from ^J3j We leave it to the reader to show that, with the 
calculated decay rate r ~ 5 x 10 _40 s _1 , even if one assumes that the entire mass of 
the galaxy, M g « 10 50 M p ; resides in hydrogen in the 3d state at an average distance 
of 30 kiloparsecs, R ~ 10 56 £ p i, the detection criterion falls orders of magnitude short 
of being satisfied. Detection of gravitons from the spontaneous emission of hydrogen 
is impossible. 

Graviton production by black hole decay is in some sense more promising. The 
Hawking temperature of black holes is T^h ~ h/kGM. To evaporate lOeV or higher 
energy gravitons requires black holes of M < 10 27 M p i or M < 10 22 g. As this is far 
less than a stellar mass, such black holes would necessarily be primordial. Although 
theoretical prejudice may be aligned against primordial black holes, observational 
constraints on this mass range of PBHs are almost nonexistent |17] and in principle 
they could make up most of the universe's dark matter. On the other hand, PBHs with 
M < 10 15 g, which would emit particles with energy > 10 8 eV, would have evaporated 
by the current age of the universe. Constraints on such PBHs are much tighter due to 
the distortion they produce on the X-ray background [17] . Let us then assume PBHs 
with 10 15 g < M < 10 22 g make up the entire mass of the galaxy. 

Detection criterion (12.61) requires to an order of magnitude 



To determine / 7 over the range of interest, note that the temperature of a 10 g black 
hole is too low to evaporate anything but massless particles: neutrinos, photons and 
gravitons. By equipartition one expects approximately equal numbers of these species 
to be radiated, but the higher spin of the graviton suppresses its emission compared 
to the neutrino and photon. Page's detailed calculations [18] have shown in this case 
that / 7 = .02. For PBHs with temperature above the electron emission threshold, f 1 
drops only to .01, so it can be considered a constant ~ .01 over the entire mass range. 

To estimate Ma/e 7 , first note that from the photoelectric equation we have a 2 k 2 = 
(uj/ujRy) — 1, where the Rydberg frequency ojr v = h/(2ma 2 ). Thus ak w (e-y) 1 ^ 2 
when the graviton energy e 7 is in Rydbergs. The mass spectrum is expected to be 
dominated by the smaller holes[T7] of M ~ 10 15 g, which emit gravitons at about 
10 7 Ry or ak ~ 10 3 ' 5 . If we assume all the holes are in this mass range, the (ak) 5 
in the denominator of the ionization cross section reduces it to only a ~ 10~ u hG. 
The high graviton energy e 7 ~ 10 7 Ry ~ 10 _20 e P 2 is similarly unfavorable in criterion 
(16. ip . If we again take M g 10 50 M p i and R ps 10 56 £ p ;, we find that the number of 
detections over the source lifetime for a Jupiter mass detector is 



N d ~ 10 



-5 




Gravitons 



12 



Properly one should average (Ma/e 7 ) over the mass spectrum of PBHs, or equiv- 
alently the lifetime of the source t s in Eq. (12.41) . but the estimate we've just done 
should be adequate. Larger holes give a larger cross section, but their lifetime r ~ M 3 
is so much enormously larger than the age of the universe, that the detections per 
age of universe is negligible. We should also mention that our ionization cross section 
is not valid in the relativistic regime, but a proper calculation may only make things 
worse; if so this result does indicate that even if primordial black holes exist, they do 
not produce a sufficient quantity of gravitons to be detected. 

A more definite astrophysical source of gravitons is gravitational bremsstrahlung 
produced by electron-ion or electron-electron collisions in stellar interiors. Gravi- 
tational bremsstrahlung has been studied by numerous authors (see Weinberg |14j 
and Gould[19j). It is calculated by much the same procedure as in electrodynam- 
ics, however there the radiation is primarily dipolar, whereas as discussed in the 
Appendix, gravitational radiation is quadrupolar. In this case one calculates ra- 
diated gravitational energy in analogy to electromagnetic quadrupole radiation, by 
taking the the square of the third time derivative of the moment of inertia tensor, 
Dij{t) = fiXi(t)xj(t), for reduced mass (i. In the zero-frequency limit, one finds that 
the energy radiated per unit frequency is: 

dE 8G 2 4 . 2 

— — = — (i v sin v, (6.3) 
aw 5n 

where 9 is the scattering angle of the electron. 

If two scatterers with charges e a and e& undergo Coulomb scattering with impact 
parameter b, the deflection angle will be 9 m e a eb/{b(iv 2 ). The maximum frequency of 
emitted radiation will be u> max ~ v/b ~ (iv 3 9/(e a eb). The important point is that the 
frequency of emitted radiation depends on 9. Integrating the above to u max yields for 
small angles 

g,- 86 ^ , (6.4) 
57r e a e b 

The gravitational power radiated per unit volume is found simply by multiplying this 
expression by the number of collisions per unit time and volume, or 

~dQ 



P = v^2n a nb J E-~rr<jil. 



a.b 



where n e and rib are the number density of scattering species. In the nonrelativistic 
regime da r /dQ can be taken to be the differential Rutherford scattering cross section 



da e\e\ 



dQ 4v 4 (i 2 sin 4 9/2' 

For a Maxwellian velocity distribution the power per unit volume works out to be 

5 v 
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Multiplying P by the volume of the star gives the power radiated in gravitational 
bremsstrahlung. We note that this formula differs somewhat from Weinberg's Eq. 
(10.4.31), which contains a lower-limit 9 cutoff, In A, due to the divergence in the 
Rutherford cross section. His expression, however, neglects the fact that the upper 
limit of integration for uj depends on 9 (and vice versa). When carried out as we 
have, the In A does not appear. 

In any case, we present this derivation primarily for illustrative purposes because 
Eq. (16.31) is only valid for low-energy electrons. Stellar interiors, where T > 10 7 K 
{kT > KeV), require the opposite extreme; and it is in this limit where the Born 
approximation used in the calculation of the ionization cross section is valid. Because 
Gould has treated in some detail the astrophysical regime, we here merely quote his 
results. He finds that for the Sun, L 7 * ~ 7.9 x 10 14 ergs s" 1 implying that / 7 * ~ 2 x 
1CT 19 . For white dwarfs he finds L lw d ~ 10 4 L 7 *, but for a typical white dwarf, L W( i ~ 
lCr 2 L*, implying that f 1W( i ~ 1CT 13 . For neutron stars he gets L ry ns ~ 10 cvcjs s 
The total thermal luminosity of neutrons stars is less well determined, because much 
of the luminosity is due to rotational energy. Recent observations [20j . however, give 
a thermal luminosity L ns ~ 10 31_34 ergs s -1 , implying f ins ~ 10~ 8 , very roughly. 

The number of graviton detections from bremsstrahlung is now easily estimated. 
If T ~ IKeV, then ak ~ 10, which gives for the ionization cross section a ~ .OlhG. 
With M and R as before, and making the assumption that all the mass of the galaxy 
resides in Main Sequence stars, Eq. (16. ip yields 

A" 7 ~ KT 5 , (6.6) 

ruling out detection of gravitons, but if a substantial fraction of the galaxy resides in 
white dwarfs or hot neutron stars (for the latter: T ~ 10 9 ; ak ~ 100), 

N iwd ~ 10; N ins ~ .1, (6.7) 

giving a faint hope. Unfortunately, the uncertainties in the various quantities preclude 
more precise answers. In this regard we point out that in this calculation and the 
previous, the numbers are not significantly changed by going to the mass and radius 
of the observable universe, due to the inverse square law; this is a manifestation of 
Olber's paradox. However, one can do better by merely parking the detector at 1 AU 
from the source, which yields 

iV 7 * ~10 3 ; iV w ~10 9 ; N jns ~ 10 7 . (6.8) 

At the tidal disruption radius r T (M^/MjY^Rj 30i?j, these numbers can be 
pushed to 

Ny* ~ 10 7 ; AT W ~10 13 ; ~ 10 11 (6.9) 

For a neutron star thermal lifetime of roughly 10 5 yearsj^D], the last figure gives 
potentially 10 6 graviton detections per year. 
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All these numbers are, however, somewhat optimistic. First of all, current es- 
timates of the number of neutron stars in the galaxy is ~ 10 9 , for a mass fraction 
~ 10~ 3 . Furthermore, detection criterion (16.11) assumed that the upper limit for the 
lifetime of a source was r = M/L. In real life, for a main sequence star r ~ 10~ 3 M/L; 
for white dwarfs r < 10 _5 M/L; and for neutron stars Tth < 10 -8 M/L. Thus, the 
white dwarf result in (16 .7p . as well as the MS and NS results of Eq. (I6.8P are also 
ruled out absolutely. At this stage, the only possibility for detecting gravitational 
bremsstrahlung appears to be putting the Jupiter-mass detector in close orbit around 
a white dwarf or neutron star; the latter might result in as many as 10~ 2 detections 
per year. 

A surprising, but perhaps significant, potential source for gravitons is photon- 
graviton oscillations due to passage of photons through the galactic magnetic field. 
This mechanism, analogous to neutrino oscillations, was first discussed by 
Gertsenshtein[2I]. If in interaction Hamiltonian (14.41) contains stress-energy due 
to a background magnetic field B, as well as to the magnetic field b of the photon, 
then the linearized wave equation for the gravitational wave h^ v contains both these 
fields in the source term. Moreover, must also satisfy Maxwell's equations in 
a gravitational field, which couples the photon magnetic field to the gravitational 
wave. Simultaneous solution of these equations leads to a wave-packet that oscillates 
between the photon and graviton states with mixing length 

L = (2/G 1/2 B), (6.10) 

independent of wavelength. (Here B is the component of the magnetic field perpen- 
dicular to the direction of propagation.) For B on the order of a Gauss, L comes 
out to be roughly a megaparsec. Furthermore, the oscillation period shows that if a 
single photon travels a distance D through the field B, it will emerge as a graviton 
with probability 

P = sin 2 (D/L) = sin 2 (G 1/2 BD/2). (6.11) 

Taking the magnetic field of the galaxy to be B m 5fiG with D the distance to the 
galactic center, yields P ~ 10~ 15 . This is equivalent to / 7 ~ 10~ 15 in the previous 
calculations. For UV photons with e 7 ~ 2 Ry the detection criterion gives 

iV 7 ~ 10 5 . (6.12) 

The situation improves if one places the detector 1 AU from a neutron star with a 
magnetic field B ~ 10 12 . We take D ~ 10km. Then P ~ 10 -14 and the detection 
criterion gives 

N y ~ 10 8 , (6.13) 

assuming e 7 ~ IKeV (which may be somewhat high for an estimate of neutron star 
surface temperature). If one put the detector at the tidal disruption radius, the 
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number of detections could be raised to iV 7 ~ 10 12 . Taking into account factor of 
1CT 8 already mentioned, we might expect 1CT 1 detections per year for the Jupiter- 
mass detector. 

In reality, however, the outlook for the Gertsenshtein mechanism is less positive. 
Large magnetic fields will result in electron-positron pair production, which produces 
an effective index of refraction that lowers the speed of light. This in turn introduces 
a coherence length above which the gravitons and photons will no longer be in phase; 
the Gertsenshtein process is quenched. Dyson has recently shown[22j that the quench 
condition is DB 2 uj < 10 43 , implying for neutron stars that u < 10 13 . This is below 
the ionization threshold for hydrogen; neither is it of great astrophysical interest 
Evidently, neutron star gravitons produced by the Gertsenshtein mechanism cannot 
be detected. 

7 Discussion 

We are now in a position to address the issue that we have so far avoided: given 
that the calculation of the cross section did not involve quantizing the field, what is 
our justification for claiming that a click in the above detector constitutes detection 
of a graviton? To see the answer, consider the upcoming generation of gravitational 
wave detectors, which are expected to be able to measure a wave amplitude h ~ 
10 -20 at frequencies of approximately 10 3 Hz. From the expression for the flux / 
above, this corresponds to roughly 10 7 eV cm~ 3 or about 10 26 eV per cubic wavelength. 
At an energy hu per graviton this amounts to approximately 10 38 gravitons per 
cubic wavelength. In other words, LIGO and its successors will never be said to 
have detected anything like a single graviton. The picture just described is one of 
a classical wave, which is what LIGO will have detected, as designed. On the other 
hand, the fluxes for the sources we have been considering amount to < 10 _18 eV per 
cubic wavelength, far smaller than the 10 eV imparted to an electron in an ordinary 
photoelectric experiment. Such a situation cannot be reconciled with a classical 
picture of a wave, because there is insufficient energy at a given location to ever eject 
an electron. Hence, if the ejected electron is registered, we are entitled to claim we 
have detected a single graviton. 

This last "if," however, is a large one. Throughout we have assumed an ideal de- 
tector, of one hundred-percent efficiency, the mass of Jupiter. This is not reasonable. 
One must therefore examine the physics of the detector itself. For an ejected electron 
to be registered by a sensor, the sensor should be located within a mean-free-path 
of the ejection point. If the mean-free-path A for an electron in this case is deter- 
mined by the ordinary Compton scattering cross section ~ 10 -24 cm 2 , the density is 

■"■We have recently learned that Zel'dovich arrived at a similar conclusion regarding the quenching 
of the Gertsenshtein mechanism in 1973. See Ya. B. Zel'dovich, Sov. Phys. JETP 38, 652-653 
(1974). 
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~ lgcm~ 3 and the temperature T ~ 10K, then A ~ lm, meaning that a fraction 
\/Rj ~ 10~ 6 of all ejected electrons reach the surface. This essentially rules out any 
of the above scenarios that were not ruled out by other factors. 

One does not wish to compress the detector, because the mean-free-path scales 
as £ 3 , but it does suggest building detectors with £ < A. As one example, imagine a 
honeycomb structure made of silicon, whose interstices are filled with hydrogen. Both 
the silicon and hydrogen would become targets (with presumably comparable cross 
sections) and the question becomes, How large can one make such a detector before 
self gravity causes structural failure? A straightforward calculation shows that the 
maximum radius of such a detector is give by 



where e is the structural strength as a fraction of the bulk modulus, v s the velocity 
of sound, and p the density of the material. Suppose p = 0.233 g cm~ 3 (i.e., 0.1 the 
density of silicon), e ~ 0.03 and v s = Ax 10 6 cm s -1 . This yields i?~4x 10 8 cm with 
a mass of M ~ 5 x 10 25 g, about .OlMEarth- With such a device the detection criterion 
cannot be met for any of the sources discussed above and detection of gravitons is 
ruled out. 

This result, however, does not absolutely exclude detection of gravitons; one can 
imagine filling the solar system and beyond with tiny detectors. At this point, though, 
the possibilities go out of sight. 

Before that point, we must address two other issues. The first is noise. Any de- 
tector needs to be shielded against background noise. Two serious noise sources are 
neutrinos and cosmic rays. The cross section for the interaction of neutrinos with 
matter is about 10 _45 cm 2 , or at least twenty orders of magnitude greater than the 
gravito-electric cross section. In a typical white dwarf, neutrino emission exceeds 
photon emission|23], meaning that 10 13-14 neutrinos are emitted for every graviton. 
Therefore, without shielding, one would expect 10 33-34 neutrino events for every gravi- 
ton event. A shield should be thicker than the mean- free-path for neutrinos, which 
for materials of ordinary density amounts to light years. Such a shield would collapse 
into a black hole. Unless one can find another way to discriminate against neutrinos, 
this appears to make detection of thermal gravitons impossible. In light of this result, 
we do not pursue shielding against cosmic rays, which would activate the detector 
material, inundating it with secondary particles. 

The second issue we have ignored brings us back to the philosophical side. We 
have assumed that a click in the detector amounts to a detection of a graviton. His- 
torically, as mentioned in ^31 Millikan's "proof of the existence of photons rested 
on measuring the slope of the graph of the photoelectrons' energy versus frequency, 
i.e., the determination of h. If one insists that the existence of gravitons is not es- 
tablished until h is fished out of the data, then further obstacles immediately present 
themselves. At the very least, one requires many more detections, in order to plot 
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gravito-electron energy vs. frequency. Unfortunately, in contrast to Millikan's situa- 
tion, we here do not have a monochromatic graviton source. To deconvolve the signal 
from the source spectrum presents additional hurdles. 

In sum, we can say that to detect a single graviton was a priori going to be a 
difficult proposition, but it was not obvious that it was fundamentally impossible. 
Although, as we stated at the outset, we have found no basic principle ruling out 
graviton detection, reasonable physics appears to do so. Perhaps the most interesting 
aspect of the investigation is that it leads to some fairly subtle physics, which, as 
the Appendix shows, has caused significant confusion in the literature. Certainly, if 
a "no graviton" law appears elusive, we do feel entitled to predict that no one will 
ever detect one in our universe. 
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Appendix: Spontaneous Emission and Hamiltoni- 
ans 

The extra factor of one-quarter in the squared matrix element (14.71) surfaces during the 
resolution of a matter concerning the spontaneous emission of gravitons by hydrogen. 
Specifically, Steven Weinberg in his standard text Gravitation and CosmologyftQ 
employs semi-classical methods to compute the spontaneous decay rate from the 
3d2 — > Is state of hydrogen by graviton emission and arrives at one answer, whereas 
the authors of the well-known Problem Book in General Relativity and Gravitation\15\ 
use a field theoretic approach for the same problem and arrive at an answer 2.1 x 10 5 
times larger. This is a significant discrepancy, particularly in that for the electromag- 
netic case, the field-theoretic and semi-classical approaches are known to give identical 
results. Although spontaneous emission did not directly figure in the body of the pa- 
per, it does provide an important check for our calculations because transition rates 
must necessarily satisfy detailed balance, i.e. conservation of energy. Resolving the 
discrepancy also leads to some fundamental issues about the Hamiltonian employed 
in computing the matrix element. 

Consider, then, the spontaneous emission of a graviton by the decay of a hydrogen 
atom from the 3d to the Is state (which one might conceivably contemplate as a gravi- 
ton source). In the semi-classical approach to conventional quantum mechanics[H], 
one does not use a Hamiltonian to compute the spontaneous emission rate. (If the 
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photon or graviton does not exist before emission, it is unclear what the interac- 
tion is.) Indeed, there is no classical analog to spontaneous photon emission, but 
in the semi-classical approach one nevertheless constructs one by imagining a cur- 
rent flowing from the upper state to the lower. One assumes the current density is 
J = pv = pp/m. Next one identifies p with the probability density e^*^ taken over 
the initial and final states, such that J = e/m^P^a- As in §|H one assumes that 
p = mv = imujv. Then J = ieuj^lr^ a) and quantization of the field is avoided. In 
classical electromagnetism, the total power radiated by a dipole is ^|-|J | 2 , where I Q 
is the total current, so from the spontaneous transition we have for the total radiated 
power: 

Ak 2 e 2 uj 2 r 

'%,A>J\r 2 . (a.i; 

If we now say that each transition emits a quantum with E = hu, where u is the 
transition frequency between the two states, we can reinterpret (lA.lj) as the transition 
rate 

r = A (A.2) 

There is no reason to believe this derivation, except that when the stimulated 
emission and absorption are independently computed, the imposition of detailed bal- 
ance leads to the Planck formula and, miraculously, a proper field theoretic calculation 
also leads to the same result. 

Thus as in §H we again copy the procedure for gravity, with one important differ- 
ence. The power emitted by a simple electric dipole is P = (2/3)d 2 , where d = ex 
is the dipole moment. By analogy, one would expect a gravitational dipole to emit 
P ~ (mx) 2 . Due to conservation of momentum, however, mx = p = for an isolated 
system and gravitational dipole radiation does not exist. The lowest radiating mo- 
ment is the quadrupole, and like an electromagnetic quadrupole the radiated power 

:.2 

goes as Q . With the substitution Q ~ ex 2 — > y/Gmx 2 , for an x that varies harmon- 
ically in time we have P ~ Gm 2 uj & x A . 

More precisely, when the source dimension is much smaller than a wavelength (the 
dipole approximation) one finds [H] that the emitted power is: 



P 



2Gu 6 m 2 r . 1 



D* j D ij --\TrD 



(A.3) 



where the moments are exactly those defined by Eq. (14. 6 j) and 2/5 times the bracketed 
quantity is exactly the average given by Eq. (14.71) . Eq. (1A.3|) is essentially the square 
of the third time derivative of the moment of inertia; the spatial average results in 
the subtraction of \\TrD\ 2 . The important point is that all the matrix elements are 



3 

now seen to be the same 
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The spontaneous emission rate is now found by plugging (jA.30 into (1A.2D . For the 
3d2 — > Is transition the normalized hydrogen wavefunctions are 

and the integrals can all be evaluated analytically in terms of elementary functions. 
Weinberg [T4] apparently calculates this transition by such a direct substitution (see 
his Eqs. 10.8.5-10.8.7) and arrives at 

023/^^3 ol3q3 ri 

Y(U2 - Is) = - - R 2 = f-±£ m 3 aV = 2.5 x lO^V 1 , (A.5) 
V ; 3 7 5 15 (137) 6 ft 2 5 15 h V ' 

where to get the second equality we have used the exact transition frequency uj = 
Ah/9ma 2 . This gives a lifetime T" 1 = 4.0 x 10 43 s. 
By direct substitution, however, we obtain 

T(3d2 -> Is) = T(3d0 -> Is) = A-?m a aV = 4.9 x lO^V 1 , (A.6) 

which is about 2.2 x 10 4 times larger than Weinberg's and roughly a factor of ten 
smaller than the Problem Book's. Apparently Weinberg's result merely contains a 
numerical error. 

We can nevertheless convince ourselves that the general procedure is correct, in 
particular Eq. ( 1A.3I) . by checking detailed balance. We assume the spontaneous 
emission rate given by Eqs. ( 1 A. 21) and ( 1A.3I) . The rate for stimulated and emission 
and absorption follow from (14.21) and (14.81) . The crucial point is that detailed balance 
will only satisfied when the factor of 1/2 in the matrix element (14.51) is replaced by 
1/4 (see Eq. flA~T2|) below). Then: 

- n m 2 /Au 4 <|A/>. (A.7) 



The intensity / for a gravitational wave is I = u 2 h 2 /8irG. Demanding that the 
spontaneous plus stimulated emission rates equal the absorption rate requires 

e -Hu,/kT {Q^^j + Gcu ^ < | A .|2 >= ^2^2 < | A .|2 > 

or 

/ H = ^^TT' (as) 

showing that detailed balance is indeed satisfied. Because the matrix elements can- 
cel out from both sides of this equation, this procedure cannot provide a check that 
they have been correctly computed but it does show that if Eq. ( 1A.3I) is the correct 
expression for spontaneous emission, then T a b s in flA.7j) and | < a\H\b > \ 2 in (14.81) 
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are the correct expression to be used in the calculation of cross sections. 

The reason for the correction can be found by examining the Problem Book's 
field-theoretical calculation. A proper field theoretical derivation is the only utterly 
convincing way to obtain the spontaneous emission rate — and the only independent 
verification of the semi-classical results. Why then does the Problem Book fail to 
produce the semi-classical answer, which one would expect it to do? The crucial step 
in the authors' calculation is to write Hamiltonian (14. 4p as 



H 



(87rG) i /2 m£!_ (A9) 

m 



Here, 0^ = h^j v G. They have also chosen the same components of the polarization 
tensor as we did in §H1 implying that = piPj/m. 

Next, in the standard field-theoretic manner, is decomposed into a series of 
harmonic oscillators by the use of creation and annihilation operators. (The 87rG is 
introduced to normalize the number operator to the energy density; see BR.) 

To get T, they compute < 3d\H\ls > for the spontaneous emission of a graviton. 
This is carried out in by the standard field-theoretic prescription with the result that 
in the dipole approximation 

f = ^i <1 '^i M> i' (A ' 10) 

(The conceptual difficulty about using an interaction Hamiltonian for spontaneous 
emission is obviated in field theory by saying that the interaction is with vacuum 
fluctuations.) 

The remainder of the solution consists of a sophisticated averaging procedure over 
all spin states via Clebsch-Gordon coefficients and the Wigner-Eckart theorem. For a 
final result the Problem Book gets T = 5.3 x 10 _39 s _1 , which as mentioned is about 
10 times larger than our result. Nevertheless, a simple field-theoretic calculation 
performed with Eq. (14.81) . does yield exactly the semi-classical answer, which shows 
that, not only Weinberg's result, but the Problem Book's must be incorrect. There 
are several minor mistakes involved, which we discuss in BR; here, however, we do 
point out a fundamental error that has been made: 

In general relativity, coordinate transformations and gauge transformations are 
the same thing, and it turns out that the Hamilonian (14.41) is not gauge invariant. In 
setting en = — e22 and eu = 621, the Problem Book authors have chosen to work in 
the so-called transverse-traceless, or TT, gauge, in which the polarization directions 
are purely spatial and orthogonal to the direction of the wave's propagation. It is 
a textbook exercise to show [IB] that to first-order in the TT gauge, two particles 
initially at rest remain at rest under the passage of a gravitational wave. In other 
words, nothing happens, which suggests that the TT gauge cannot easily be used 
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for this problem. Properly, physics should be calculated in a locally inertial frame 
(LIF), in which the effects of gravity are absent. Otherwise, in the TT gauge for 
example, not only are "non-inertial" forces present, but the laws of electromagnetism 
are significantly modified as well. In a LIF, such effects are second-order, typically 
smaller than in the TT gauge by a factor of v 2 , where v is the velocity of particles in 
the system. 

In an LIF, the energy density p dominates over the momentum and pressure, 
and so the only non-negligible component of the stress-energy tensor is Too = p. A 
straightforward coordinate transformation from the TT gauge to the LIF (see BR) 
yields for the Lagrangian density 

Clif = -hj'pheJTxW, (A.ll) 
Thus, the LIF interaction Hamiltonian is 

H L i F = - j Cd 3 x = -u^whe^x^. (A.12) 

We see though, in the case pi = imuxi, the locally inertial Hamiltonian itself has 
exactly the same form as the TT-gauge Hamiltonian in (1A.9I) . but that each term is 
smaller by a factor of ~ 2. Thus, it is not entirely surprising that a TT-gauge calcu- 
lation leads to almost the locally inertial result. Nevertheless, the LIF Hamiltonian 
( 1A.12I) is the easiest one to use in the computation of the cross section. One can carry 
out the calculation in the TT gauge (see BR) but in that case one needs to include 
the electromagnetic stresses in the stress-energy tensor. 
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